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A procedure was developed to aid in the testing of a 
two-sided composite null hypothesis about the mean of a 
normally distributed random variable for situations where 
either the population variance is known or unknown. The 
procedure was designed to eliminate the requirement for 
iterative type solution techniques normally used in deter- 
mining the acceptance or rejection region of the subject 
hypothesis. This thesis provides guidelines, curves, and 
tables which will aid in testing a two-sided composite null 
hypothesis. Provisions were also incorporated into the 
procedure to permit testing of hypotheses about the relative 
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. Specified lower (d.) and upper (d 


SYMBOLOGY 


Probability of a Type I error. 
Probability of a Type II error. 


Critical value defining the rejection or acceptance 
region centered around Wy (o known case) or zero 

(o unknown case). 

The lower bound (C,) and upper bound (C on the 


) 
critical region. R 


5 
xX. 
i=1 + 


n 


The true population expected value (mean). 
A specified value of the true population mean. 


R) limits for 
the true population mean (u). 


Population variance (o = population standard 
deviation). 


Sample variance used to estimate o%, specifically; 


xX. - raya 


iy 
n-l 


a n ( 
6o. 2 ) 
1i=1 


sample size. 
Cl yn 


Oo 


6 Yn 


io} 
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ee : , Oo 
Coefficient of variation ; 
n (aT? 
A specified value of the CV. 


The null hypothesis of interest. 
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H = The alternative hypothesis to Hy 


P = The tail probability that a value equal to or 
more extreme than the observed X would have 
resulted if Hy were in fact true. 


xr = Desired maximum allowable eo. of the true 
mean from some specified value (nu in terms of a 
portion of the population 5 ae deviation 

§ (i.e., |[u-ull < ro). 

a6 = The statistic 


Yn(x—- u,) 


‘Fa 


DF . = Degrees of freedom. 


r = The non-centrality parameter of the non-central 
Pdistripution, A’ = x7n. 


6 [x] Normal cumulative distribution function. 
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I. INTRODUCTION 


A. BACKGROUND AND PURPOSE 

A procedure was developed to aid in testing a two-sided 
composite null hypothesis about the mean of a normally dis- 
tributed random variable. The intent was to eliminate the 
need for iterative solution techniques when solving for 
critical values manually. In situations where computer 
assistance is inappropriate or unavailable, there is a tendency 
for the analyst to shift from the interval test requirement 
to the computationally easier two-sided simple hypothesis 
(i.e., that the mean is equal to the midrange value of the 
interval). By employing the procedure proposed in this thesis, 
the time required to test a two-sided composite null hypothesis 


by hand can be significantly reduced. 


B. SCOPE 
Procedures were developed for the following test requirements. 
lL. wFor Testing, That qe So tis d, 
For situations where the population variance is known, 
a method was developed for testing that the population mean 


is between two values (dy, and d,)- 


23° ‘For’ Testing That ens Stee 


For situations where the population variance is unknown, 
a method was developed to test that the absolute value of the 


deviation of the population mean from some specified value (HO) 


a a d Ps are 
ae2 nd bias oF Bagoleved’ ‘eed “ exabe90r 
7, Ay ree 


Yo meson ecty gueda i mozoqul £ tie 


aN a 
= > > ‘ ‘ 

< - RRS , a2 
aul er o 
Je Le; “ . 

— ° “ 
> it ¢ 
a r 

‘ ii o of 


ny 620OG i Pee Sv wb wsayv ce 3 


~~ een 2» 2a07 poisaeT 39 
ef 28S 


Re a a a 


it$siveog exit ened aliescneh nisi aN) 
edy faye) ie on +32. begalavad (ein Boe 
igh bam 15) anetav x 
= 7 M5 


Noss aLi@og 
s 
ith od ane as 


sere 4 #073 oui 


is less than or equal to a specified portion (r) of the 
population standard deviation. By a slight modification, 
this procedure can be used to test if the coefficient of 
variation is between zero and some specified value. 

Each procedure consists of using relevant situation 
parameters (i.e., 56, n, 6, a, etc.) to determine the critical 
region using eieher interpolation curves or tabulated data. 
The procedure was designed to be versatile with respect to 
using Type I (a) and Type II (8) error level control, or 


using tail probabilities (P,) for drawing statistical 


T 
inference. 

To simplify use of the procedure, the paper is 
organized so that the main body contains only 1) minor 
theoretical context, 2) detailed procedural guidelines, and 
3) example applications. The detailed derivation of the pro- 
cedure is provided in Appendix C for readers so inclined. 
Appendices A and B provide curves and tables used to find the 
associated rejection or acceptance regions and the tail area 
probability. The curves are provided to simplify critical 
value interpolation for what is considered to be the majority 
of applicable parameter values. However, the total breadth 
of parameter ranges considered in this paper is covered in 
the tabulated data of Appendix B. It should be noted that 
all explanatory figures in the main body of this paper are 


schematic in nature and do not represent scale drawings. 


10 


td snolaitteoe eds Y Hr of bean 
ata, cdelateial th — ‘baa o2se ded 8 

noliapa Le anawod ers pain 0 “esetence ‘exbevera fk 

tesitsias ad ealmetoh 02 088 ) a) a 9 so 

pees bodaludss +8 aevioe aqksafouzsant ‘aedghe | 

o3 sodgees ftiw sliseaiiy od og banpieed asw 

.. 30 ,fortnon Level wT Va) Tt oat bas (o) I 3 

igotjeivese palveip see Ue) asizitidedong tal 


et zeqeq edt ,etebeooag afd Bo sar yiitente c 
uontm (tf y¥ieo enipgoon “bed nian sais jen? oe, 
Pris eonilebiog { g3wbaserg. Ballegveb (s able etetindated 
ais Yo 4 soigsvizeb beltisved emt enoiss21 lege 
Spaifoni on. axebtis, BOT > 2tbweageh as bebivialll 
ie Olt oF heen zelda? Dan Sava siivesy € Sas A 
fens 158) ef2, Ses enolges sphetgesos 3p) roisaetex 


ei 


isnpigizo yIiiamle od Sebiverg exe aerare ont - aE 
(217i ee, sat HR..P2, hazvebleden, 4b gsiw, 208 sto 
Joosd Latos eft ,vevewon , eautay satameTEg ¢ 
er evon ai tedeqg eid? of becV5obiados 2 

(ats be9on, od blued ay 28 pani 
eis 78 98q ig te ybot niem ‘odd nd seaupld 
,Bpad ~- elspa teenengex Jon ab baa MS 


6 
9 Oe, hm 


C. SYMBOL AND HYPOTHESIS TESTING CONVENTIONS PECULIAR 
TO THIS THESIS 


Whenever possible, symbol and concept conventions which 
conform to standard hypothesis testing terminology were 
employed. However, because of requirements which are peculiar 
to this procedure, it was convenient in some instances to 
use terminology that may be new to the reader. The following 
symbol and concept conventions apply throughout this Thesis: 

1. The Null Hypotheses (H) will always be bounded 
intervals such as Pa = es <-0705 j= wel < ro. All alterna- 
tive hypotheses are the complements of such intervals. 

2. The level of significance is customarily a pre- 
assigned upper bound on the size of the critical region 
(i.e., that subset of the sample space reserved for rejection 
of Ho): This will hereafter be referred to as a error level 
control. | 

3. By reversing the role played by the null and alterna- 
tive hypotheses, it is possible to define an acceptance 
region for the originally formulated interval null hypothesis 
with a specified preassigned upper bound (8) on the size of 
a type II error. This procedure is defined as 8 error level 


control. 
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II. RESULTS AND DISCUSSION 


A. CASE OF KNOWN POPULATION VARIANCE 
1. General Remarks 

In order to test the null hypothesis that the popu- 

lation mean falls between two values (i.e., d 


< < 
1,2 UE ap), 


an iterative solution of equation 1 is required. 


e 
i] 
hk 
oO 
+ 
+ 
‘oO 


where, 


o is the known population standard deviation, 


C is the value of the statistic x which defines 
the critical region centered around 


_ 4, + 4: 
Ug 2 ’ 
a is the probability of a Type I error, 
Sate 
6 is equal to a a and 


6 is the integral of the standard normal probability 


density. 

To determine a rejection region for the subject 
hypothesis, an initial value of C would be guessed and the 
right side of equation 1 calculated. A new C would then be 
used iteratively until the desired accuracy in the value of 


a had been obtained. This solution process is not difficult 


or time consuming with the assistance of a computer. However, 


when performed manually, such a process may require more effort 


than is warranted by the scope of the project. 
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The procedure developed in this thesis will provide 
the individual with a method for testing a two-sided com- 
posite null hypothesis which circumvents the previously des- 
cribed iterative solution technique. The procedure derivation 


is provided in Appendix C. In general, it consists of 


6¥n 


normalizing equation 1 by reparameterizing the terms ae 
CYn 
ej 


" 


and as the parameters "a" and C*, respectively. Appro- 


priate values of C* are then obtained according to the pre- 
scribed procedure which follows. 
2. Procedure Description 
The procedure for testing that dq, ss d, is out- 
lined in Figure 1 (for ao and 8 error control) and Figure 2 


(for estimating P This procedure was developed for 


.): 


Situations where symmetric critical regions are of interest. 
A summary of the procedure follows. 


a. For Controlling The a Error Level 


First calculate a = a If the value of "a" 


is less than or equal to 1.7, the associated value of C* 


is obtained from Figure A-l of Appendix A. The lower (C_) 


L 
and upper (Cp) bounds of the rejection region are then 


defined using equation 2. 


*O * 
G2) GSS! nex . fy Chere ier Seo where Uy Sythe & 


L ie) Ja R e) /a 


If the calculated value of "a" is greater than 
1.7, the relationship between "a" and C* becomes linear and 
equation 3 is used to define the rejection region (Ky is 


obtained from Figure 1). If x < C, or x > C, then reject 
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the hypothesis with a probability equal to q of rejecting 


when in fact the hypothesis is true. 


iS (a + K,)9 (a + K_)o 


b. For Controlling The 8 Error Level 


First calculate a = ney If "a" is less than or 


equal to 3, the associated C* value is obtained from Figure 
Aaa) sae S196) (Or iFigure A=ZeCha6 < a'< 3)- Of Appendix A. 


The lower (C_) and upper (Cp) bounds of the acceptance 


L 
region are then defined using equation 2. 

If "a" is greater than 3, the relationship 
between a and C* becomes linear and equation 4 is used to 
define the 8 error level acceptance region. If Ch < x < Cp, 
the hypothesis that the mean is between dy and dp is accepted 
with a probability of falsely accepting equal to 8. 


(4) Co =ae pincats Se and CL =u_ t+ 


L fe) Sa R fe) ee 


(a - K_)o (a - KJ9 


c. For Estimating Tail Probability (Pr) 


The. exact value of the’tail probability (P,,) can-~ 


T 
not be calculated since it depends on the actual value of the 
population mean (u). However, the procedure depicted in 


Figure 2 will provide an estimate of the maximum possible 


value of P,, for a corresponding value of x equal to or more 


ak 
extreme than that observed. The procedure for estimating Pr 
consists of first calculating the values a = sill and 
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ie} 


Go to Append.A Go to Append.A(Fig.2 
Fig.1. for .2¢a<1.6 or Fig.3 
for 1.6<a<3). 


Calculate C at desired error 
level curve and intersection 
of vertical line projected 


Calculate C & Cp: 


Tega Cea 


upward from calculated value ‘a: pee ee 
a - 
“Uy* yn - 


and C 


RTS 
CS > Cig For 8 control 
R Yo Yn E 


If X<C or X>Cp then 


reject H(d)< u < dp) 


then accept 


R 
H(d) su < dp) with an g 


probability of falsely 
accepting 


with an a probability 
of falsely rejecting 


APPROPRIATE Ke VALUES SYMBOLS 


eve 
d, & d.=Lower and upper test interval] 
values. 
6 = -5(dp- qd, ) 
Va -5(d) + dp) 
C,= Defines lower bound of critical 
region. 
C,= Defines upper bound of critical 
region. 


Figure 1: Algorithm For Determining a or 8 Critical Regions 
for Known o 
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Yes Are No 
2.8 and Ch € 5. 
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Yes 


* 
aC 23.52 


Go to Append. A 
(Fig. 4 for O<as1.5 
or Fig. 5 for 
].5«a¢2.3 ) and deter- 
mine point where 
vertical line from 
'a' value and 
horizental line 
through C* intersect. 


Calculate Pr as 
P_=1- eLatC*]+ofa-C*], 


where © is the integral 
of the standard normal 
density. 


Calculate Pr from: 


Sie hast. 
Po efa-C ] 


Interpolate between 

Seb to determine 
(Note: these curves 

14 plotted log 

and should be inter- 

polated accordingly). 


Figure 2: Algorithm For Determining The Tail Probability 
For o Known 
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- These values are then used to estimate 


ti 


Pr by employing the interpolation curves of Appendix A or 


by direct calculation using the equations provided in Figure 2. 
3. Range of Parameter Values Considered 
The procedure itself is suitable for all possible 
parameter values. However, a restriction on the range of 
parameters covered by the curves and tabulated data was 
necessary. The range of parameter values was selected ina 
manner which would provide adequate coverage for the majority 
of conceivable real-world situations. 
a. Curves 
1. For a Error Level Control. Parameter ranges 
covered are a(.01, .025, .05, .075, .1, .15, .2, and .25) 
and a(0 to 1.7). 
2. For g Error Level Control. Parameter ranges 
covered are 6(.01, .025, .05, .075, «1, .15, .2, and .25) 
aad a(.2to 3). 
3. For Estimating Pp. Parameter ranges covered 
aera (0 to 2.8) and C*(.04 to 5.2). 
b. Tables 
1. For a Error Level Control. Parameter ranges 
covereduvare- oi.001,... Od, (3025), 055," 5075;, Gali, «LS, .2, and 


-25) and a(0 to 3). 
2. For §8 Error Level Control. Parameter ranges 
sovered are @(.001;,...01, 2-025, .05, .075, .1, .15, .2, and 


.25) and a(0 to 5). 
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4. Example Applications 
a. Example 1 

In the production of a certain structural element, 
a new size has been established as marketable. The expected 
size of the new element is to be between 16.5 and 17 cm in 
order to satisfy engineering and management requirements. 

For example, if the expected size is below 16.5 cm the 
rejection rate of unacceptable elements will be too high. 
If the expected size exceeds 17 cm, more material is used 
than is necessary. 

For the first 25 units produced, the average 
size is 17.03 cm. It is known from previous experience that 
the element size is normally distributed with a standard 
deviation equal to .2. In addition, no functional relation- 
ship between the mean and standard deviation of the element 
size has been observed in the range of 10 to 20 cm. 

The decision maker is interested in determining 
if the hypothesis that the expected size is between 16.5 and 
17 cm can be rejected with a Type I error equal to .05. The 
test is donaucted with 6 equal to .25, Uy equal to 16.75, and 


(225) (5) 
2 


a equal to or 6.25. Since a is greater than 1.7, 


the rejection region is defined using equation 3, specifically: 


C, = 16.75 - S:-22 Ft Sho) t2) = 16.434 


and 
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CR = 16.75 + .3158 = 17.066. 


Since x is neither greater than C, or less than C,, the 


L’ 
hypothesis cannot be rejected with a Type I error equal to 
OD. 

b. Example 2 

From a sample size of 81 shipping containers 

selected at random, it is necessary to know if the expected 
container weight (which appears to be normally distributed) 
is between 75 and 75.5 kg. The sample average weight and 
variance are 75.22 and 1.44, respectively. A Type II error 
equal to .01 is desired. Since the sample size is relatively 
large, a reasonable estimate of the population variance can 


be obtained using the sample variance. The test is conducted 


with 6 equal to .25, Uy equal to 75.25, and a equal to 


(25) (9) 
ee 


a C* value equal to .0726 is obtained from Figure A-3 of 


or 1.875. Since the value of a is less than 3, 


Appendix A. Using this value of C*, the values of Ch and 


CR obtained from equation 4 are 

C= Low 2on= {.07283) (0-2) = 7.25 6: 00975 = 75.24 
and 

Cp SF Ta 2o OOF = 7526 


Since x is not between Cy and Cpr the hypothesis that the mean 
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is between 75 and 75.5 cannot be accepted with a Type II 
error equal to .01. However, it may be possible to accept 
the hypothesis at a somewhat higher risk (8) of falsely 


accepting. Trying 8 = .05, the new C* obtained from Figure 


A-3 of Appendix A is .35 and accordingly Ch = 75.2 and CR = 75.33% 


Since x is between Ch and Cpr the hypothesis could be accepted 
with a probability equal to .05 of falsely accepting. 
c. Example 3 

For the situation of Example 1, the decision 
maker wants to know for what values of a the hypothesis 
Geos iu < L7 could be rejected. ‘The answer can be found by 
estimating the maximum tail probability associated with the 
observed x using the_ procedure described in figure 2 with 


(x- uo) Yn 


a = 6.25 and C* = = 7. Since a + C* > 3.5, 


the estimate of Po would be Pr = $[.25-7] = .227. Therefore, 
an a probability greater or equal to .23 would be required 


before the hypothesis could be rejected. 


B. CASE OF UNKNOWN POPULATION VARIZANCE 

1. General Remarks 

A procedure for situations of unknown population 

variance which would parallel that of the previously dis- 
cussed variance known procedure cannot be developed without 
several stage sampling such as Stein's Procedure (reference 1). 
However, the previously developed methods can be extended 
directly to the case of unknown variance provided the null 


hypothesis is scaled by the standard deviation (see below) 
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and the noncentral t distribution is used in place of the 
normal. Details follow. 


a. A procedure was developed to test the hypothesis 


eee 


Bie Ue | eSex ol Or equivalently < r where 


uw is the true population mean, 


U is some specified value of interest near 
which the mean is to be tested, and 


r is the portion of the population o for which 
the test is to determine if the true mean 
iS within ro units. of Hos 

b. The procedure can be used to conduct tests about 
the displacement of the coefficient of variation (CV = Tuy? 
from zero. 

2. Procedure Description 

The procedure for testing that |u- Ug! sire ds 
provided in’ figure 3. For testing that CV < CVO, figure 4 
should be used. Finally, the method for estimating Pa is 
provided in figure 5. 

As the sample size gets large, the ratio ns converges 
to the test value of r. To take advantage of ie eee. 
the curves for a error level control are presented in terms 
of the value Ee instead of just C. This provides an indica- 
tion of the atte of r and the sample size n for which the 
6? unknown curves become close to those of the 0? known case. 


However, Since the method of estimating P,, requires repeated 


T 
table entry and interpolation, the critical values presented 
in the tables of Appendix B are in terms of C directly. A 


summary of the procedure follows. 
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a. For Testing That |u- ug ee) 
Select the desired a or 8 error level. Consult 
the corresponding curves of Appendix A or tables of Appendix 
B. If the curves are used, multiply the value obtained for 


oe by the sample size (n). If the tables are used, determine 


n 
the corresponding value of C directly for the subject values 
of r and DF. Calculate the statistic | 
(x- UW) vn 


eee 


If the test is to control the a error level then reject the 


hypothesis, with a Type I error equal to a, if the value of 
r is greater than the C value determined from the a curves 
or tables. For a test to control the g error level accept 
the hypothesis with a Type II error equal to 8 if the value 
of x is less than or equal to the C value obtained from the 
1-8 tables. 
b. For Testing That CV < CV, 

To test the null hypothesis that the coefficient 
of variation does not exceed CV, at the a level of significance, 
one need only apply a role reversal technique to the above 


procedure. The set equivalences 


illustrate the fact that the hypothesis of interest is one 
of the alternate hypothesis treated previously, and can be 


Managed by reversing the role and interpretation of a and 
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8 error level control from the previous test for 
ju-ull < ro. 

1. For Controlling The a Error Level. Obtain 
the C value corresponding to the subject values of r (equal 
to a) and DF from the table of Appendix B associated with 
the valie i-a. Calculate the statistic 


xv/n 


AZ 
Oo 


If r < C,--xreject..the hypothesis that Cv < CV, with a Type I 


error equal to a. 

2. For Controlling The 8 Error Level. Obtain 
the C value corresponding to the test values of r and DF 
from the curves of Appendix A associated with the value 8 
Wee. , ror 6 = .0L, consult the curves for an a =..01). 
Calculate the statistic r. If f > C, accept the hypothesis 
with a Type II error = 8. 


c. For Estimating Pr 


The actual value of Por associated with the observed 


value of xr will depend upon the true value of 
= dea) 


Oo 


To obtain an estimate of the range of possible values for Pre 


the procedure detailed in figure 5 is used to calculate a 


value for P,, corresponding to r set at zero and r set at the 


ay 


test value (r). The actual value of Pr will be between these 


two extremes. 
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Calculate 


For aerror control 


Go to the appropriate figure 

of Append. A for the associated 

a level(Fig. 6-15) or cables 

of Append. B for parameter 

values not covered by curves 

of Append. A. Project vertically 
from the r scale to the appro- 
priate OF curve then horizcntally 
to the C//n_ scale. Multiply by 
vn to obtain C. 


If r> C, reject the »ypothesis 
Ho with a probability = a of 


rejecting when in fact it is true 


Figure 3: 


Algorithm For Testing The Hypothesis 
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For g error control 


Go to the appropriate table 
of Append. B associated with 
the desired 8 level(i.e. for 

B =.1,go to the C value table 
associated with a ,P2,1-8 =.9). 


Project down the corresponding 
r value column and across the 
respective DF row to determine 
the correct C value (Note: If 
the r value being tested jis 
between two of the table's 
columns ,interpolation is 
necessary.) 


If r < C accept the hypothesis 
He with a probability = 8 of 


: accepting when in fact it is 


false. 
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Calculate r= 


xn. 


For controlling a error level 


Go to the appropriate 8 table 
of Append. B (Tables 3 —» 41) 
which correspondsto ang value 
equal to the desired a level 
Cie.0 FOr a=, go" to. the. C 
critical value table associated 
with 1-8 =.9). Project down 
the corresponding r value 
column and across the respec- 
tive DF row to determine the 
correct C value. [ Note : If 
the r value being tested is 
between two of the yable's 
columns , interpolate for C by: 


(Test "Left 


High""Low? 


If r < C, reject the hypothesis 
CV «CV, with probability «2 


of rejecting when in fact it is 
true. 


Figure 4: 
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For controlling 8 error level 


Go to the appropriate a figure 
of Append. A (Fig.6——15) which 
corresponds to ana value equal 
to the desired 8 value(i.e.for 
8 =.01, go to the a=.01 curves 
of figure 6 or 7). Project 
vertically from the r scale to 
the appropriate DF curve then 
horizontally to the C/vn 
scale. Multiply by Yn to get C 
[ Note : For parameter values 
not covered by the curves ,go 
to the tables of Append. B 
which provide C directly , 
instead of C/vn.] 


Ifr>c » accept the nypo- 
thesis CV<CV with 


probability 8 of accepting 


when in fact it is false. 


Algorithm For Testing The Hypothesis CV < CV, 
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Calculate r 


First calculate the estimate of P. for r set equal to zero as follows: 


Go to Tables 3447.Pick a P-(i.e..5) and go to the associated table 


Project across the respective DF row and down the appropriate 
r column to obtain the correct C value. (Note: If the r being 
tested falls between two columns,interpolation is necessary. 


Go to a lower Po value 
table(72e.°<2), A 


as that table No 
been checked on a 


previous step? 


Yes 


The P- for this r is between the current table Po value and the 


previous table. For a closer approximation interpolate as follows: 


A 


Now repeat the 
process to get 
an estimate of 
Pr for r= ir. 


Have the 


estimates of Pr fo No 


both r=0 and r=r bee 
alculated? 


Yes 


The true P. lies between the two estimates for r=0 and r=r 


Figure 5: Algorithm For Estimating P,, For o Unknown 
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3. Range of Parameter Values Considered 


Two general levels of r values, along with selected 
incrementation, were chosen to parallel the range considered 
applicable in the majority of situations. For example, in 
most cases where the test of interest is that lu-ugl Seo 
yr would range between zero and 1. Conversely, most tests 
concerning the coefficient of variation would require that 
-143 < CV, < 1 (or equivalently, 7 > xr > 1). However, the 
eritical values (C) for the total range of r (i.e., © to 7) 
provided in the curves and tables can be used for either 
test. The specific parameter ranges covered were: 

a. Curves 

Curves are provided for a error level control 
for both low and high values of r. 

1. fox, Low x Values (0 < © < 1). Curves are 
provided for a equal to:.01, .025, .05, .1, and .2. 

2). HOD High is Valuesi.(2 <5. < 7)« 4 Cuxves are 
provided: for -a, equal to,..01, .025, .05, .1, and .2. 

b. Tabulated Data 

The same tables are used to present critical 
values for both a and 8 error level control and for estima- 
ting Poe To obtain C values for 8 error control, the tables 
associated with a and P,, equal to 1-8 are used. The values 


T 


covered by the tables are a, P or 1-8 equal to .0Ol, .025, 


Yd 
MO Serie cle: weet Qaptetods: partis Aine si DipintrerOlsio cl pe hsiOtn Meco ae OO pa ueoLD puand 99°. 
For each value of a, Pre and 1-8, three tables are provided 


to cover the range of r from 0 to 7. 


27 


7) o ey ae sa 
~ 4 mz Se, , : og bs 


eiqnexe 203" Books se siz big 


pe 1 Bibow noize irev te 
a ’ 
re lore 6YV 
i ay we i- t + tu 
werd e WD Fine 
\ =) « ; - ot -}; val BOs & 
’ 
; t t bs $3055 
. mf 4 
: 9 g — > et. 
a 
: ’ [ %O7 , pS cows) | 
" 4 ; 
' oat ed 
_ 297, 8V ipld: (bas we bh 
: 


brit «tl. 4«f0s- 4 Sn. Loy oF, ‘Taupo 5 308 


— 
_ i ae 
wer) (yt + > £ {oti eB > Aout t0¥ x ) 


a al \ 
. Bas’ ts (200 ISG. VIO ee 4 Incps. cM, 
atad | bagasnaat” 
Lgodsiue’ dneteaxg GF seu era: aoldas onuba ett . 
pid r . i , ie , 4 [Py Py _s 


vpelses t0ot fas LowineD water note 4 be, 2 ason 


et wit ans od aren 


aN 
vy i. a ae 


<5 


wy 


4. Example Applications 
a. To Test That Juul < x16 


It is necessary to determine if the output of a 
new system is between 11.75 and 12 units. The standard 
deviation in the output of the new system is believed (with 
confidence) to be between 2/3 and 3/4 of that of the old 
system (o = 1.5). Using the conservative estimate of the 


improvement in the system deviation (i.e., 3/4), the test 


C120. =8 1 14.7 5!) 
C2) C875) (15) 


For the first 15 units, the values of the sample mean and 


parameter r would be set equal to or 2.1L 
deviation were 11.68 and .3, respectively. The desired level 
of significance (a) is .05. 

The test is conducted with Uy equal to 11.875 and 
r equal to 2.517. Using figure A-10 of Appendix A, 7 is 
determined to be .6, and accordingly, the critical ae 
equals (.6) (W15) or 2.324. Since r is greater than C, the 
hypothesis is rejected with a Type I error equal to .05. 

b. To Test That CV < CV, 

A relationship is known to exist between the mean 
and standard deviation of the purity of a certain chemical 
product. A new method of processing is being considered to 
improve the expected purity of the product. The decision maker 
is interested in determining if the same degree of quality 
control would be obtained with the new process as with the 
current method which provide a standard deviation equal to 
-18 of the mean (i.e., a coefficient of variation equal to 


-18). 


28 


7 7 are. 2 


uA Ae +a ane ita 
oe he: pure yom 
a %60 YugrHo add FL onteas308 od Rann, 
branes ant. ddttey 80 Bos cet id) ‘ae 
djiw} Seyallfed ez shill coll Wen ond te Seest0 ucts 1a 
bio edd 20 salts ao ' BNE bas e\s needed" od oi 
ois 2 oJamides svistwvrsenoo and ecieu” «tes re 
$asy érid eee «+8, ft) nolseiveb nedeve, edd as | 


vehy too eae oF taupe tok Pt © bivow tee | 


bis chem siqase efx 46° neulav’ ens eg the ai seen 
level Serleesd edt  .ylevigoegess: ,f5 BRE aa irk . ve 
<Goe O22 ta) is . 

itt re. fd rr Laup ys Ay ow ee at deod nee | ‘ay 
A xibueqggh 30° OL-4 exonkh SAlag  “<Siae 
euisv Lavizizo ond , ¢benibiécoe baw , a. bdo t } 
eid. nets nésnexz éi x asnle -bLE,8. 20 (ZEN 
20. of Inups voitzxa ft egy? s dire SesSeted ez ak 
5 2 v3 35nT tat or, 

isem eit Aeewsondt J8ixne od awond el gidenasdeted a7 . 
[xo tuiedo nittuss a fo yoebyog’ ele 86 te lseiven = aug 

oo hexehlones poled #¢f on tsdecotd to hodzer wer Poe 
toxan ftoteisel adt southony py Fo ysixwy Lereonae nasa 
Jilara Se sexvseb emsa ads 53 paLebarsedab: “nb 

eta do iw Be ssepotg wan ott da be bends ait 


te 


- 
aire 


ae ee oe 


os Lenpeliteaives bsebaaga ‘s obivoug. 4 dw 'b on . 
ot a oe to banat is yp Aid 30 | 
‘ i i p ee ) 

Ms a ee A im. at 


4a 
’ 


Since a relationship between the mean and standard 
deviation is known to exist, and only 15 batches of the pro- 
duct have been produced by the new process, the population 
variance is considered unknown. The mean purity of the first 
15 batches was 4.6 parts per million with a sample variance 
equal to .25. The decision maker is concerned that if he 
does not convert to the new method, he will be rejecting, with 
a probability (a) equal to .05, an opportunity for a better 
product. The rejection region associated with a value of r 


equal to 5.56 (i.e., ) is found by consulting the critical 


se 
-18 
value table with 1-8 equal .95. Interpolating between the 
table columns for r equal to 5 and 6, a value of C equal to 
17.69 is obtained. The value of r is calculated to be 35.63. 
Since r is greater than C, the hypothesis that the new process 
provides a quality control equal to or greater than the 


current method cannot be rejected with a Type I error equal 


ie, +05. 
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APPENDIX A 


CURVES FOR DETERMINING «OR 8 CRITICAL VALUES AND TAIL PROBABILITY 
(P 
T 


I. GENERAL 


The curves provided within were developed to simplify 
procedural calculations in lieu of the tabulated data of 
Appendix B. They do not cover the total parameter ranges of 
the data in Appendix B but were designed to maximize interpo- 
lation precision while covering what is felt to be the majority 
of applicable situations. For those cases which are subject 
to parameter values outside the range of the curves, Appendix 
B should be used. Figures A-1 through A-5 are relevant to 

2 


the o” known procedure. Figures A-6 through A-15 deal with 


the ae unknown situation. 
nD ACCURACY 


The accuracy achievable in using the curves vs. using the 
associated tables depends on the interpolation prowess of 
the individual user. However, to provide a general idea of 
the accuracy which may be expected, the following statements 


are considered germane. 


A. ae UNKNOWN PROCEDURE CURVES 


a. C* Values For a or 8 Error Level Control 
Accurate interpolation to the third digit is possible 
ee ei,. C 


equals C* -005). 


* + 
true interpolated ~ 
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b. Po Values 


Because of the semi log scale and interpolation 
between curves, the accuracy in the interpolated Pr values 


is expected to be ~ + .02 for .1 < Pr <))-Oppand +)).005. (tor 


Pr od OF Pr SE giy 


B. 57 UNKNOWN PROCEDURE CURVES 


a. Curves For The Range of r From 0 to l 
Interpolation accurate to the third digit is 
possible. 
b. Curves For The Range of r From 1 to 7 
Interpolation accurate to the second digit is 


possible. 
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C* Values For a Error Level Control With o Known; 


Figure A-1: 
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C* Values For 8 Error Level Control With o Known; 


Figure A-2 
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C* Values For 8 Error Level Control With o Known; 
(6.45 var <3) 


Figure A-3 


34 


ed (aes 7 


ee ee eee 


~ 


£ 


6.5 Rl Bet teh 


se% 


Dice U5 


Figure A-4: Tail Probability for o Known with 


35 


r r 


1 


; 
t -* 


% 


P 
} 


g, 


a2iw nwoet.o 10%, yo 


ak 


2 


ae 


R« 


Figure A-5: 


—o0o4 
na— 


oh Oo NWOPON 


nm 
on 


Tail Probability For o Known With l. 


36 


Das Ales, 2.8 


vet 


CVn 


0) 


sey Raia Cakes ipa Sa ame 


r 
Figure A-6: C//n Values For o Unknown With a, Bey Of 1-8, =) 01; 
(G siaths 


Oo 
er 
ry, 
w 
- 


37 


+07) eae 


: tal ue = vt: ~ - —t - : ; | SAGs 
LH ated fae ye" i 5 a ate 


-O1; 


Pa or 1-8 = 


(ee) 
38 


C/Yn Values for o Unknown With a, 


Figure A-7 


CERO SARee eRe OD Dees 
08 SESE Bee eee 


rrr 
ees oe! O07 40 Bee 
A 4 A 
Pity) Pi At SOC c Oe oe 
SRREE BERS’. Bee eee 


Sane (0SGS GGGS6 CRGRGGeGGe GReRe) Geen Coeee Genes Sees ene 
EECCA 


= SHITE Hien 


See ices tects cased pened cumee Ceeaeey 


a 


rrirr) 
ae ee See eee 
BERR Os PERE eee Pee eee 7 ane aa 
tI SERS Be OR es Ree” eens Pee ee T 
Co SS A 5 SD 


ee | 
aS BERS eee” ae es Bee 
r +p——— Jane f 
SRUP= 28808 SESE Bee eee [ 
Bl SRLS A HH ek 


Cee 
3s =—— 20888 SRR ES Pee PRRs Sees Cee eee Ae 
eS 


N 
H 


MY 


LANE 
VAAN 


| 
\N 


i 
H-A 


\\ 


\\ 


\ 


rs SS NOL 
TS YH SO 
SRS OSS es eee 
DEPRU SEES eae eee eee 
es oe | 


Figure A-8: C/VYn Values For 
(xr 


Unknown With a, P 


, or 1-8 = .025; 
1) q 


[A Q 7% 


a4 
mae Di OE ee, 
» md Oe & Be 
ni ear a 
< > ae 2 
oie eases 
; Seer r siete 
a >. io oo 7 hesdnbaplpg ote hogar pe ob teak Aen yo 
4 aoa i? ° i ery Sas ~ ined oe 
s : Pe Fbe 6408) 6 y2 whinge o hin, 
; = rf tit l lonedyn tees 
S or MM pane wid be S55 
; 7 = a See Se 
> ' ; : i eel wee eeie a4 Seaeerort pers 
i +3 . 8 2 SoS ee Ba aoe hae Te be: 
t JRE. } 
i , } } 
$e paeeeudae Sn 
bt 
’ > ~~. 
‘ i 
‘ = po 
| 4 t- — ies 
- . ; ; : $o44 peed oo pte 
° ’ ataey tagoos 
‘ 
se - 
¢ -. 
- . a 
. f 4 . “ 
i 4 
- a ; 
» ow nme 
— rs 
{ * i 
j 
Peeves eee ee ee eee ee ee ee ee ere 
a es Hy oth ire oeress verre) Ste 
> 4 n oom bee om Gem Ay to St ter f 
¥ + aie feed ar angen 


6 ree ew 


a a. ee 


is ~ (i 2 
i an ’ ; i 
0. = Rel to gad yo A22W mwonigy: a 
i eek 
vy) 7 ¥ _ 


ad 


, 
pes 


7 
ry, 
i 
, 


20 


18 


16 


14 


C/vn 


= fF HS ra oa 2oenn oS rs 
OA cEnanar ora} me 
— 7 = fA + tt + 
ip ses 7 + 
is é 
‘ a" sep am t 
Ae ae ia 
4 : Sree eee re = 
Z Pet 1a55 CoE TH Tas r 
3 is +H f I poet am 
35 yas nam sas 0SaS8 
A y i rt r m1 och tr t 
‘a saanaee osenee _oaneae +I 
z 1 soe qo opt rr see =} = t 
A tot ttt - - - > = 
t t t r ra sqoncaac ceo r 
cam 2 as = I 
T T Tttrt + : 
t +t t + ret 
+ ae 
Br = 
a0 SOSS8 BSR HHH + T 
+t T 
ose 


] Z 3 4 5 6 7 
r 


Figure A-9: C/yYn Values For o Unknown With a, Pr jer 178 = 
(ie iT) 


40 


pit 


inv 


ee 4 


re 


ont 
2 
oC 


nw 


t i" 


2 


A 


\F 


ie te ; 
= a = ' f : yf ; 5 oy 
re ; aS ~~ i ia ay a : 5 
a > - 7 fo 68a 'S 


3.0 


2.8 


2.6 


2.4 


2x2 


2.0 
C/yn 


1.8 


1.6 


1.4 


Aw 


Figure A-10: 


C/Yn Values For o Unknown With a, P 
(a? 1) 


41 


7? 


or 1-8 


=e 


ae 


teat Peet 
acoteNates! 
Ley: 333 eats 


. 
ape wi hereap— p pantie 


t ¢ pie a 


. 
Wha pee hae se ana gre eet | ; Bisaie toes 
ies Senn oe Sere & SS era Bae. SO 

3 tel Sarhet teas Seiteteeeaees 


Yt rer te coe bape, 
baee sey ++ J 4 tae (4044845 
op ee ers - 


eb 4 pee were T > 
+ 


20 


18 


16 


14 


C//n 


Figure A-1l: 


Fae ca Axa 60 
SS EAS AIA ee 
= i ay so zeae —o co 
; pe jp sdaeer 25 py ast soe - Seanet 
=. ‘ae ss FAN roa 
ce = s tt tt 
‘as ral ci peasaneaes : { 
zZ be srs eases 
2 + = +H 
aes Zz isoee nea t + con 
Ty — 2 i> 31 —_ 
im > asaune wae iS 
eee attics Hee 
= : 
= : z 


2 3 4 
if 


C/Yn Values For o Unknown With a, P,, or 1-8 


(eg) 


42 


-05; 


ae td. 


ah 


3.0 fesse sere es egeeeeg==ieeepesi se et fee 
Oe eee ge lesiay ga pany uersay ee 
ge = = oo eae eS =e 


8 ee eee ee Eee Ye. E | 
26 fa fA 


C/ vn 


Seer LLL ILE 
12 age f ZE 
Eee LLL” N 


Seg 


Figure A-12: C//n Values For o Unknown With a, Pen eae = a1; 
(eos) 1) 


43 


7 - = >| 
Liat Chas hig el 


hea 


pS Sl SEAS 


a 


fey 


ae Tr) 


‘ 
4 


— 
} 
‘ 
' 
i 
' 
7 ‘a 
* 
| 
— pe 


a 


Y 


=3.. 
ha 


JDF=2 


22 


sam 


Bs eens 


20 


18 


16 


14 


12 


san chee 


C//n 


7 


6 


4 


3 


2 


r 
Values For o Unknown W 


Figure A-13 


Eh) ay, Pre or 1-8 = .1l 


a 


C/ 


Gaps) 


44 


18) >Re 
See ee 


EE AA LAA 

ALIS 

aac EEA Ae 

ELLE AAA 

- AA i 
a OMY 
= naa yyy 
Jas yy Oe A 


= eee 


“ 


r 
Figure A-14: C/vYn Values For o 
< 1) 


at bow 


; = ~ ht 72 


ay) 


a 
= 
fea p 


C/vYn Values For o Unknown With a, 


(2 ys) 
46 


Figure A-15 


) “AatW mortal o 
Xo 
ac 


ae : 


e . 


TABLES OF a OR 8 CRITICAL VALUES AND P 


APPENDIX B 


Table B-1: 


a 
a 0,001 0,01 0,025 0,05 0,975 0,1 0,15 0,2 
a 

0,0 3529056 2,57583 2,2414 1,95996 1,78046 1,64485 1,43953 1.28155 
0.1 3,3967G6 2,588F 2,25254 1,96972 1,78934 1,€5306 1,44677 1, 2879f 
OleSis 3,326.34) 1260426) 2.26629 1, 9.8282 1.800371, 66327 1.45569 2, 29596 
0,2 3,35258 2,62564 2,28519 1,99°55 1,81567 1,67748 1,%682 19017 915 
Ceizom eos says E52O. 24 FOR93 2.01971 “Les St1 169550 eee 2? 2.32152 
0,3 SM ANS) AVES SIWR 20. IS711 2,.04/S05 1905851 1, 70785 1,50367 1.33903 
0.35 3,46146 2,71871 2.36928 2,074296 1,88565 1,74992 1,5276€47 1,35964 
0,4 3,59459 2.75732 2,40496 2,10699 1,919625 1,771.73 1,55249 1,3F33 

0.45 3,5%981 2,79867 2,44364 2.14286 1,95093 1,803°2 1,58159 1,40991 
0.5 3,59654 2,84223 2,4R48A4 2.10148 1,9°666 1,83875 1,6136 1,43937 
O.55 3,64434 2,88751 2.52811 2,22245 2,9258 1.876208 1,6483 1,47153 
0,6 3.69:2228 293409 2.57304 2,26539 2,0671 1,916€12 2,68547 1,58623 
0,65. 3.74193 2,98166 2,61927 2.30996 2,11024 1,95795 1.72485 1,54%327 
0.7 Bia? Dit3. .3' 02996. 2.66652 2.395589. 2.4:5494 2.00047 1, 76619 1,5.8244 
0.75 3,84091 3,097877 2.71453 2,40275 2,20082 2,04642 1,89092 1,62352 
0,8 3 890/66 v3 6127916) 24 76392 2.45047 (22477312, 092'52 1,852366 1, C6626 
O85) 3.9405 (3,17741 2,8227112 2,4988 2,29542 2,13955 1,89933 1, 71044 
,0.9 3,9904 3,2270% 2,86141 2.5476 2,34371 2,18733 1,94598 1,755a4 
9.95 4.04034 3.27579 2.91093 2.59675 2,39247 2,235F8 1,99347 1,89224 
1,0 4,09031 3,32663 2.9606 2,64615 2,44158 2,28447 2,0415 1,8%947 
71,05 4,14028 3,37653 3.01038 2,69573 2,49094% 2,33359 2,099N6 1,89736 
1,1 4.19027 s.42E6€46 3,06024 2.74544 2,54949 2,38297 2,13961 1,94577 
1.15 4,24626 3.47642 3,11014 2.795294 2,59018 2,43257 2,1f70974 1,9946 

abia?s 4,29026 3,52639 3.16007 2,84511 2,63997 2,482921 2,23769 2, 94372 
1.25 4,34026 3,57637 2,21903 2,89502 2,689°2 2,532 2 2enary 24 0932 

1.3 4.39025 3,62636 3,260N1 2,94496 2,73972 2,5°185 2.337903 2,14965 
1.35 4,44025 3,67636 3.39999 2.99492 2,78965 2,63175 2,3P°6°3 2,17233 
1,4 4,49025 3,72635 3,35998 3,0449 2,83961 2,6f1F8 2,4367 2,2471 
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APPENDIX C 


PROCEDURE DERIVATION 
I. KNOWN POPULATION VARIANCE 


A. ERROR LEVEL CONTROL 
The a rejection region associated with the hypothesis 


that the mean is between dt and d, is defined by 


ao = P[x < eri + PL x CL] 


where the rejection region is chosen to be symmetric around 


Uy (see figure C-l). Since the suprema will occur when 


yu equals qs or das the critical region for the specified upper 


bound on the Type I error (a) can be written as 


Gite Pisce < C, |given ie = Uowl + P[x > Cp | given es ugrs). 
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Since the sample mean (x) will be normally distributed, 


this can be expressed in the form of equation (C-1l). 


{C)-(u- 6)} Co-(u,- 
L Mo" 6) ym ae re R bee 


Beta e phy 2 < 


where z = Lee (in 5)1 ~ N(0,1) 


Substituting Cc = Le and Cp = Hot, equation C-1l becomes 


(C-2) a =P[ Ze KE {s-C)yr ] + P[ zZ> ( SON | ] 
or 
auzeep 2a = “Ay 44 - of of7 + Cyn 
oO (oJ Oo oO 
By reparameterizing, the number of parameters necessary 


for solving equation C-2 can be reduced from five to three. 


Letting a = os and C* = cy equation C-2 becomes 
* * 
(c-3) a= 1+ afa-C] - ofatC] 


Using the result that Cr. = Wy7 or CR = Ho tc and 


* 
C= eee equation C-4 is obtained and is used to define the 
n 


rejection region. The hypothesis is then rejected, with a 
probability of falsely rejecting equal to a, if the sample 


mean is less than Cr or greater than Cpe 


(C-4) CAS) uaa 
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For each value of a there is a value of "a" such that 
the relationship between "a" and C* in equation C-8 is nearly 
icinear (i.e.; C* = K, +a). Bor a > 1.7 the relationship is 
approximately linear (starting with 5 significant digits 
for the worst case with a = .25, and improving as "a" in- 
creases) with slope one for all values of a < .25. The value 
of Kp can be obtained by considering that if @{a+c*] ~ ll, 
then equation C-3 simplifies to equation C-5. This can be 
rewritten to provide Kp = o° (1Saie When this linear relation- 


ship occurs, C. and Cp can be obtained using equation C-6. 


L 


Mesa) =e elas Cpe 1 or o= tla = cc ] 


B. 8 ERROR LEVEL CONTROL 
In order to determine an acceptance region for the 

hypothesis qd, ORES dp, with a probability no greater than 

8 of falsely accepting, the alternative hypothesis pe Se > 6 

is tested. The a rejection region of the alternative hypothe- 

Sis will be the 8 acceptance region of the null hypothesis 

(i.e., d, < u < d,) and will be defined by C, < 2c 


L 
(see figure C-2). 


R 


The suprema will occur where u = dy. or d,- Thus, the 


upper bound on the Type II error (8) for the hypothesis 


d, <u < dp, can be written as 
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sta < C,| given that u = p_- 6] 
ueH yo ae 
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This can be expressed in the form of equation C-7. 


C i - 
(C-7) «@ (H, ) = P[z A R ae _ ppze (C) soy 


[x -(u, -8) ya 


fof 


where z= 


By substituting C_ = roe and Cp = Tete and noting 


L 
that z will be normally distributed (i.e.N(0,1)), equation 


C-7 reduces to equations C-8. 


(c-8) a(H,) = af eae BMEe Sef ee 
or 
a(H.) = a+C]-efa-C] 


Equation C-8 can be reformulated to resemble equation C-3 
which will permit use of the same solution technique as was 
used for a error control. Specifically, by multiplying 


equation C-8 by -l, adding 1 to both sides, and setting 
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1-o (H) equal to a', equation C-8 reduces to equation C-3. 
The solution of this equation with a' = 1-a (H,) will also 
be the solution for a' = 1-8(H,)- 


For each value of 8 there is a value of "a" such that 

the relationship between "a" and C* in equation C-8 is nearly 
linear (i.e., C* = Kota) - For a > 3 the relationship is 
approximately linear (starting with 3 significant digits 

for the worst case of 8 = .01, and improving as a increases) 
with slope one for all values of 6 > .01. The value for the 


constant (Kp) can be determined by noting that when 


O[a+cCc*] ~ 1, equation C-8 reduces to equation C-9. 


* * -] 
(G=9) a@ = 1-6 = ofa-C] or a-C =90 fl -8 ] 


ae h = = 
and C =a - Ke where Ke =o [ 


The value of Ke for 8 error control will be the negative of 
the value of Kp for a error control. To allow the same Kp 
value to be used for both a and 8 error control, equation C-10 


would be used for 8 error control when this linear relation- 


ship holds (i.e., a > 3 for this procedure). 


| (a= Ko )vn (a= K.) /n 
(C-10) Ch = Wy Sere and Cr Ss Ho + eg 
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C. ESTIMATING TAIL PROBABILITY 


The actual value of Pr for an observed value of x will 


depend upon the true location of uy within the interval d, 
to d, (see figure C-3). Since uw is unknown, the exact value 


of Pn cannot be calculated. By setting u = He > Oy the maximum 


possible value of P can be obtained from equation C-3 by 
(X- uo) va 


for a, and tee ge, for C*. For large 


TE 


substituting Pr 


values of a+cC*, this reduces to equation C-ll. 


(C-11) Pr =: o{[a-c*] for atCc* > 3.5 
ae HO Ugt(tg-X) 
: P. with ae 
4 d u as 
| L 0 uytlu,-x) 
P. with u = ge 6 = dy 


Figure C-3: Dependence of Pe on Location of True u Value 
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II. UNKNOWN POPULATION VARIANCE 


A TO PEST THAR |u| < xo 
1. a Error Level Control 


The statistic 


is used to perform this test and has the density of the 
absolute value of a non-central T, as depicted schematically 
in figure C-4. The associated a error level critical region 
is defined when the non-central location is at r, resulting 
in a non-centrality parameter (\) equal to rvn. The corres- 
ponding critical value (C) is obtained from equation C-12 
and the hypothesis is rejected, with a Type I error equal 


to a, if r is greater than C. 


0 r=Vvn C 


Figure C-4: Density of the Non-central T(Absolute Value) 
Distributed # 


(Czl2 yo = PL ors, C/ given x.= ryan] 


2h 8 Error Level Control 


The suprema for the value of 8 will occur when 


A = rvn. Equation C-13 is used to determine the value of C 
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which will define the acceptance region with a 8 probability 


of falsely accepting. 
NG) iP Level 2) < Cia 


By a slight reformulation, the same solution technique 
employed in solving for a critical values can be used here. 
Substituting a' = 1-8 into equation C-12 and solving for C 
will provide the critical value weaeet ated with the desired 
level of a Type II error = 8. This equivalence is shown in 


equation C-14. 


(C-14) g = PE {Ti <c] = 1- PE |Ti >c] or PL|T| >C] =1-8°=a 


Because of the equivalence between the C values for 


a and for 1-8, the critical values for a, P and 8 can be 


Tv! 
presented in the same table provided that for a desired 
value of a Type II error equal to 8, the table corresponding 


to a = 1-8 is used. 


B. TO TEST THAT CV < CV, 
1. 8 Error Level Control 
To formulate the procedure for testing CV < CVO. 
first consider testing the hypothesis _ <r or equivalently 


CV 
Le 
| = =| S 2. SBy Setting Wy equal to zero, the null hypothesis 


can be rewritten as || < r. Solving equation C-12 with 


u. = 0 and rt = )xvn 


ig | will provide the rejection region with a 


V6? 


101 


Wedel 
ae 


i 7} a a } hy ee 
Vier ee 


oe 
ee 
a peg % 


erpindos? solsutor sige eng .nobsalumzoles sheite’ 
_etec beau ed asp eauvlaw fs92sizo 6 16% priviog, ‘ 
> tot paivicw Bas Sf+9 noissupe osnd Ant "opel 
Kexiesb ott dyiw Sevstoosan @uley’ (esieizo obit , 


nt owode 2) sonslavieps eid? +a * se73e8 ne 4 


wo ael= (Je (TER -a0 [O<cl Tl 39 ef aos: ee 


wiev 2D edy neswied socelaviupe Bay 30 eee or 
rt 


}@ 
¥ 


rN 


x4 nee. of he ..90 .o wOd Beutsy imoka a>, welt as 
bertaeb s 102 sade? bebh/owg sidsz) oie art? 4 

eres 

stinotbesyon eljat.en? ,o OF taupe 2o238 2h.e0 \ 2 


to 

Ree >. VO hehe 102 sipbeocodiq eff a3 
ylanetev viii 3° + % sleedsoqer ons panies 
cinadsoqyl Eis wid, o269 oF? sepa 9 | 
noiw nholgaups patvlo® a 


2 dote a moitoetes aa? oh <a 


Type I error equal to a. However, rejecting the hypothesis 
that a < xr, when r is greater than C, is equivalent to 
rejecting that CV > 2. This in turn is equivalent to 
accepting that CV < CV. with a Type II error equal to the 
value of a. As a result, it is possible to test the hypothe- 
sis that the coefficient of variation is less than some speci- 
fied value (CV) by using the previous test procedure for 
testing that eae < ro, provided the following differences 
are observed: 
a. Obtain the critical value (C) from the a tables 
(or curves) associated with a value of a equal to 
the desired value of a Type II error (8) and 
with r equal to a 
° 
Dien eee ris greater than or equal to this value of C, 
accept the hypothesis CV < CV, with a Type II 
error = 8. 


2. a Error Level Control 


The same general formulation philosophy applies for 


__@ error control of the test CV < CV, as for § error control. 


Again, first consider accepting the hypothesis ay <.F with a 


Type II error = 8. This is equivalent to accepting that 


CV > 2 (or CV,)- However, accepting the hypothesis CV > 2 
with Type II error = 8 is equivalent to rejecting that 
CV < ‘ (or CV.) with a Type I error equal to the value 68. 


Thus, the desired test that CV < CV, can be conducted for 


controlling the Type I error level by using the test procedure 
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for [te < ro and observing the following differences: 
a. Obtain the critical value (C) from the 8 tables 
(i.e., the table associated with the value 1-8), 
where the value of Bis set equal to the desired 


value of a Type I error (a) and r is set equal 


i 
to a. 
CVG 


b. If £ is less than this value of C, reject the 
hypothesis that CV < CV, with a Type I error 


equal to a. 


B. ESTIMATING TAIL PROBABILITY (P,) 


T 
The true value of Pr will depend on the actual value of 

=U 

=| which, under the null hypothesis, may vary from zero 


to r (see figure C-5). The range of possible values for P 


ab 


can be estimated by calculating the Po associated with the 
observed value of r with i} set at the two extremes zero and 


rv/n in the solution of equation C-12. 


P(A at O)= WH) 


an 
WOON, 
=0 r=rvn r 


Figure C-5: Dependence of Pr on the True Non-centrality Location 


103 


ae aes pee | 


ovtev ‘isolat irs 


~*~ é 
« ~ 
pig ioumen wites 
orn f 
& Co oo). ee 
orig PT eqyr 
- 
ta 
J, ? = 
ka 
} oe 
r 
4 h 
’ a | 
= 
P +88 
f 
“a - 4a? @ 
Pa 
Pa 
pant 
Ps 
“ 
al 
- 
eee 
eel Rh 


oe ED 


Ny ee ' 

DUITAMT TS. 

a 
Py... 


= , 
euleyv spay 


APPENDIX D 


SOLUTION TECHNIQUES AND COMPUTER PROGRAMS 
I. GENERAL 


The same general techniques were used in solving for the 
critical values in both the known and unknown variance 
situations. This consisted of a combined proximity interval 
search and Newton-Raphson technique. A master program was 
used to monitor and control solutions for C values at various 
selected combinations of the respective parameter values. 
Solution accuracy was controlled by stopping the iterative 
procedure when the desired accuracy in the a or 8 error 
level was obtained. The specified accuracy for the case of 
known o? was .000001. In the case of unknown o%, the 
specified accuracy was .00001 for all parameter values 


except the following: 


1. a= .01 with 125 degrees of freedom and r = .9 

2. a= .1with 90 degrees of freedom and r = .9 

3. a= .3 with 125 degrees of freedom and r = l 

4. a= .3 with 60 degrees of freedom and r = 5 and 6 
5. a= .4 with 50 degrees of freedom and r = 2.5 


The specified accuracy for these latter parameter values 


was .00005. 
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Il. a? KNOWN SITUATION 


The programs used for this case were SIGMAKWN, CINCR, 
CINCRNEWRAPS, NEWTONRAPS, and NORMPROB. SIGMAKWN was the 
Master problem which stepped solutions along the selected 
aor g values for each specified man value. Controls for 
switching solution subprograms were necessary for parameter 
values whose respective C* were too distant to permit solution 
solely by the Newton-Rapshon method. Equation D-1 was used 


during proximity interval searching (CINCR and CINCRNEWRAPS). 
(D-1) ' F(a,a,C*) = 1- @[a + C*] + Of[a - C¥*] 


Equation D-2 was employed when the Newton-Rapshon method was 


used. 


* * 
2.50663 {1-o[atC. ,]+e[a-C. ,]-a} 
PpEoVIC e Coat 7 cite pp RS ae ee ; 
oF / aC UF Exp[-.5(a+C,_1)°]+Exp[-.5(a-C,_,)°] 


Tit. 9? UNKNOWN SITUATION 


The programs used for solving for critical values when 
o was unknown were MHUSIGRATIO, MHUSIGRATIOSOLENR, CINCR2, 
CINCRNEWRAP2, NCT, IB, and TRUNERR. MHUSIGRATIO was the 
Master program stepping solutions along the selected r values 
for each specified DF. This permitted sole Newton-Raphson 


solution once the C for the first r value was determined. 
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To solve for critical values pesoerseed with testing the 
o? unknown two-sided composite null hypothesis, the equation 
P(r > C] = a (where r has a non-central T distribution) was 
solved for values of r between 0 and 1 using the method 
discussed in reference 2. This method requires use of 


equation D-3. 


xe 542 


3) ge 2) ? 
(D P[|T|>c] 7a acm 


where x =C + ( DF +C 


‘To solve the infinite series of the product of a Poisson and 


Beta distribution, the length to which the series must be 


summed in order to maintain a truncation error less than a 
specified amount (ec) was determined using the conservative 
estimate provided by equation D-4. Equation D-4 is solved 
by iterating m (the length of the series) until the left hand 


Side is less than e« (this task is performed by TRUNERR). 


(0) nat 
(D-4) E a! <6 


This method of solving the non-central T is quite adequate 


for small degrees of freedom (DF < 30) and non-centrality 
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values (A < 3) but becomes lethargic for higher values. As 
a result, a second program was written in FORTRAN which used 


the non=central T subroutine of the IMSL packages. 
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V ARRAY*SIGMAKWN sJUNK 37 3LOOPS ;ALPHO 3CSTAR ;DIVERG 

[1] ACLOBRAL VARTABLES=ALPHA(C VECTOR OF DFSTRFED ALPHA LryrL); 

C2] APASTAR (VECTOR OF DESIRED LITTLE A VALUTS);O0P7TION(O FCR 

C3] PALPHA CO"TROL,1 FOR BETA CONTROL) ;ERROR(DISIRED ACCURACY) 

C4] JUNK +(T+1),(LO0OPS+<10), (CSTAR+((C1+pALPHA )x(pACTAR ))10) 

[5] CSTAR+((pASTAP),(1+pALPI'A))p CSTAR 

{6] L1i:A+ASTARCT] 

[7] +L6x1T=1 

fs] 2L4%12>|(A-ASTAPC I-11) 

C97 L6:C+A CINCPPITWRAP ALPHA(L1] 

C10] Lu:s+«2 

[11] CSPT/ARCI311+ASTARCI] 

(12] L7+«12 

£13] £211:ALP+ALPHAL S-1] 

C 4a] NEWTONRAPS 

[15] +L7x1DIVERG=0 

C16] C+A CINCRMFRFRAP ALP 

[47] L7:CSTAPCI371<C 

C1e] t511%1(2+pALPPA )2J+JS+1 

239) C+CSTAPCI 32] 

[20] +h1x1(1+pASTAR )e#I+I+1 

[21] 37-L2x,O0PTIONe1 

[22] ALPHO+0,1-ALPHA 

[23] +L3 

f2u] L2:ALPHO+0,ALPHA 

[25] L3:ARRAY+,CSTAR 

[26] ARPRAY+ALPHO,ARRAY 

[27] APRAY<(I,(J-1))p ARRAY 

{28] "AxVECTOR | +------------------ ALPRA/C VECTOR (S )--------------- 
Vv 


V C+A CINCRIEWRAP ALP ;DELC LL 3SIGN 3F JUNK sLOOPF 
C1] sUNe<(DI7C+1), (SIGN+1), (LOOPF+0),(C+A), (LL+25 ) 
C2] +L2x1ALPso,5 
C3] C+(A-3) 
fu] £1:C+C+DPLC 
[5]. +55x1LL=LOOPP+LOOPP+1 
C6] L2:P+1+((NORMPROB(A-C))-NORMPROR(A+C))+(71)*ALP 


zal +L3x 1 LOOPF20 

[eq +L3x1\F>0 

[9] STGN+71 

[10] £3:+51x10<SIGNxF 

C11] +24 .DFLC<1 

[12] +1 1LOOPF=0 

[13] C+C-DFLC 

C1u] +21,DFPLC+OFLC+#2 

C15] Lu:NEWTONRADS 

f16] +0 : 

C17] £5 :'#*e*eDIVERGENT AT ':L03'LOOPS,FINAL C='3C3'! AND DPLC=!'3PCLC 
Vv 
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{1] 
[2] 
[3] 
C4] 
[5] 
[6] 
(7] 
Ce] 
C9] 
19] 
oheokal 


fa) 
{2] 
{3] 
C4] 
(5] 
[6] 
7] 
Cal 
{9] 
[10] 
fag 
fad 
39 
C14] 
(15] 
(16] 
[17] 


G2] 
[2] 
(3] 
C4] 
(5] 


VY NEWTONPAPS sJUNF 3~LL3;LOOP ;00nN sUP3CTEMP 
TUN Ke (LL+-12),(L00P+0), (NI VERG+9 ) 
L1:UP+#1+((~1)*xMORMPROB(A+C))+(NORMPROB(A-C))+(( 1)xALP) 
DOWN+(*((70,5)x((A4C)*2))) +0 *(0 70.5 )x(C4-C) *2))) 
23% iDAWN=z0 : 
DIVERG+1 
>0 
L3:0+C+2,59662875xUP£DORN 
>L2*1LL2LO0P+LCOP+1 
C+A CINCR ALP 
+0 
L2:+L1ix1.ERROR<|UP 
Vv 


V C#A CINCR ALP3;MELC3LL;SIGU3F;JUNK ;LOOPF 
JINE+(C+#A), (LOOPF+0), (SIGN+1), (LLI+50) 
>L2*x1ALP<0,5 
" C#(A-3) 
L1:C+C+DFELC 
+,5*1L=LOOPF+LOOPF+1 
L2:F+1+((NORMPROB(A-C))-NORMPROB(A+C))+(°1)xALP 
+L3x 1 LOOPF#0 
+L3*\F>0 
SIGN<« 1 
L3:751*10<SIGNxP 
+L4Ux 1 ERROR>|F 
>01*1L00PF=0 
C+C-DELC 
>L1,0FLC+DFLC+2 
Qus'"FIRSTC = '3C3"' WITR DPLC=' sDELC;' AYD LOOPS=';LOOPF 
>0 
LS: '*xex*DIVERGENT AT '3LL3'LOOPS,fINAL C='3C3' AND DELC=';DEFLE 
V 


V PeTORMPROB X32 
+3x1~a/(}X¥* 6,576,5l%)e 0 6,5 
+0, P+0,549,5xxX 
Z~ 0,31938153 ~0,356563782 1,7814779397 “1,821255978 1,330274429 
P<(* 0,5*¥xX)3(02)*0,5 
P+0,54+(*X)x0,5-Px( (4140, 2316419% [X)o, ayS)+, x2 
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VY C+A CINCRNEWRAP2 ALP3;LL;SIGN;3F ;LOOPF;JUNK ;DELC 
fi] AGLOBAL INPUT;C(IWNITIALK GUESS) ,DF 
{2] JUNK+( SIGN+1),(LOOPF+0 ),(LL+30),(C+A-DEFLC) ,(DELC+1) 
[3] +L1%1\((ALP20,05 )+(DF25))21 
C4] JUNK+(LL+60),(C+A-DELC),(DELC+10) 
C5] £1:C+C+DELC 
(6] +L5*% 1 LL=LOOPF+LOOPF+i 
C7] CI+LAM NCTC 
C8] F+CI-ALP 
C9] +L3x*\LOOPF 20 
(10] -+23x.F>0 
(11] SIGN 1 
C12] £23:%01*x10<SIGNxF 
C13] 2*24x.\DFLC<0,51 
[14] L2:C+C-DELC 
(14] +L1,DFLC+DFLC+2 
-[16] L4:DERVIAT+(CI-( LAM NCT(C-1E-8)))#1F'8 
(17] +L2x.DFRVIAT=0 
{1e] +L 2x \DELC<(P*DERVIAT) 
(19) WEWTONRAPS2 
[20] +0 
C21] L5:'e*x*xDIVERGENT AT '3;2L;'LOOPS IN CINCRNFWRAP2(T='s3I3'J='373!° 
{22] C+A CINCR2 ALP 
Vv 


VY NEWTONRAPS2 3LL3;CTEMP;JUNK 
{1] AGLOBAL INPUT:ERROR(DFSIRED ACCURACY) 
[2] JUNK+(LL+12),(LOOP+0),(DIVERC+0) 
{3] F+CI-ALP 
C4] L1:FTEMP+F 
{5] C+C-FtDFRVIAT 
(6) CI+LAM NCTC 
E73 F+CI-ALP 
C8] +0x 1 ERROR>|F 
(9] +L2* 1 LOOP>1 
(10] DERVIAT+(CI-(LAM NCT(C-1F 8)))#12 8 
C11] £2:*21* 1 LL2LOOP+LOOP+1 
[12] 'xxeeeNEWTONRAPS DID NOT CONVERGE(IF='323',7='303'), SWITCRED 
{12] DIVERG+1 
C14] +0Ox1KILLe1 
(15] C+CSTARCI; ds 1] CINCR2 ALP 


| 
| V PTUPPFR+LAM NC™ C3PNSSO;JUNK 3IBZFROsLL3I 
| £2] AGLOBAL INPUT:PF,CYCLE (NUMBER OF LOOPS PFOUIRFD TO 
| f2] a CONTROL TRUNCATION ERROR ) 
r3] AGLOBAL FNPUT+P=DF42sDF+NS-1 
| JUNK+(POSSO+(1t4*l SAMY), (IBZFRO+P [8 0,5), CxKeD Pet? 7+ C%2)) 
(5] TUHK«( PTUPPER+POSSOXIBZFRO), (I+1) 
[6] +Ox CYCLE =0 
[7] Li: Lh+POSSOxLA"+I 
[8] PTUPPER+PTUPPFER+(P IB(I+0,5))x*L2 
9] POSSO+LL 
(10) +61%1.(1+CYCLE)4I+T+1 
Vv 
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(11 

[2] 

3] 

C4] 

[5] 

[6] 

(7) 

Cal 

[9] 

C10] 
f11] 
(12) 
(13] 
(14] 
(15] 
[16] 
C17] 
C18] 
[19] 
[20] 
C21] 
[22] 
(23] 
C24] 
253) 
[26] 
[27] 
[28] 
(29] 
(30] 
ote) 
[32] 
33] 
[34] 
(35] 
(36] 
(37] 
(38] 
[39] 
[40] 
[41] 
C42] 
C43] 
fou] 
C45] 
C467 
C47] 
(4a) 
[49] 
[50] 


V ARRAY*MHUSIGRATIO 3A 3NS3PF 5P3C3LAMV; CYCLE 5 LAM: LAMV;LAMSQ23;XX;LAM 


ACLOBAL VARIABLES=LI7TR(VECTOR "FST \yRU/SIG{s LITTR) 
ADFVECT( VECTOR OF DEGREE OF FREEDOM) ;ALP=LOS;E2ROR=DESIRED PRPCISION 
JUNK+(I+1),CA*LITTRO11),(KILL+0) 

DFVFCT+,DFVECT 

CSTAR+((1+pLIPTR )x(pDFVECT))\0 

CSTAR+((pDFVEC™) ,(1+pLITIR))pCSTAR 
PRECISION+1+ERROR 
Li:DF+DFVFECT(I] ; 

'DATA FOR DF=";DF;' WITH ALPHA='3;ALP 

' LOOP TIME LITTR Cc 

P+DF +42 

LAMV+LITTRx(DP+1)*0,5 

CYCVFCT+LAMV TRUMNERR ERROR 

CYCLE+CYCVECT{£1] 

TIMEO*+I21 

LAM*+LAMSO2T1] 

2L5x%1\A>0 

A+0,0001 
L5:C+A CINCRNEWRAP2 ALP 

+L6*x\DIVERG=0 

C+A CINCR2 ALP 
L63:J+2 

CSTAP(LI311+DFPVECT(T] 
L11:+L13x ,\CSTAR(I;J1= 1 

2L7x\J=2 

TIMEO+I21 

CYCLE*+CYCVFCT( J -11 

LAM*+LAMSO2(J-11 

CI+LAM NCT C 

DERVIAT*+(CI-(LA”" NCT(C-1EF &)))4#1" 8 

NEWTONRAPS2 

+L 7X \DIVERG=0 

KILLI+1 

C+CSTARLI 37-1] CIMCRNEWRAP2 ALP 

KILL+0 
L7:CSTARLI;71+C 

TIME T+( (121 )-TIMEO )+60 

' ';LOOP;' 's( TIMES ;! 'sE0TFRET<-11;'° 's¢ 
2L13x%,\J=2 : 
MULT+1 

PLUSxX ASC (LITPPCL T -1150,4)4+(( 7-1 ap LITT) ) 
MULT#(LITURG TI -LITTRE TS -1 1) eC LOTTI CS -1] -L7T TPC J -27) 
115 :C+C+MULTx (C-CSTARL I 37-11) 
[13 3:+L11%1(2+pLh777TR )eJeJ+1 


+L1%\(1+pDFPVFC™)#I+I+1 
PPMPLITTR+O ,LITTR 
APRAY+,CSTAR 
ARRAY+TEMPLI?TR ARRAY 
APRAY*(I, (1-1) )pARRAY 
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V APPAY+MHUSIGRATIOSOLENR OLDCAPRSASNS3973P3C3CYCLUSLAMV sLAMEN2 SXXsL/M 
fi] AGLOBAL VAPIABLES=LTT™OR(VFOMD mrem IMI /STE|< LITT?) 
[2] ADPVECT(VECTOR OF DIGREF OF FREFDO');AlLPHLOSsP2ROR=DFPEIPFD PPECISION 
(3) AOLDCAP?(2xpDFVEC™ APRAY OF PPEVINUS C'S WITH ROW1 FOUAT, 70 
C4) a 0,0LD LITTR 1-2 AMD \-1,00L1 FQUAL 70 O,DFVFC™) 
[5] ACFSTIMATE( VECTOR OF FSTIMATED C FOR PIRET™ ELTMENT OF LITA Til BACK DF 
(6) CESTIMATF+LITTR(1) CESTIMATOR OLDCARR 
{7] DFVECT+,DFVEC™ 
C8] JUNK*(T+1) ,(A+CFSTIMATEL11), (KILL+0) 
9] CSTAR+((1+pLITTR)x(pDFVECT))10 
[10] CSTAR+((pDFVECT),(1+pLITTR ))pCSTAR 
(11] PRFCISION+1+ERROR 
(12] £1:DP+DFVECT(T] 
[13] 'DATA FOP DF='3DF;' WITH AXPHA=';ALP 
(14] ' LOOP TIME LITTR Cc! 
(15] P+DF#2 
(16] LAMV+LI™TAx(DF+1)*0,5 
(177 CYCVFCT+LAMV TRUNERR ERROR 
(18] CYCLE+CYCVECT{(1] 
(19] PIMEO+r21 
(20) DLAM+LAMSO2(1] 
[21] C#CFSTIMATF(T] 
(22] CI+LAM NCTC 
(23] DFERVIAS+(CI-(LAM NCT(C-12°8)))#1E°8 
[24] NEWPTONRAPS2 
[25] L6:J+2 
[26] CSTAR(I31]+OLDCAPR(T+133) 
[27] £11:>L13x1CSTAPL I s7]= "1 
{28] +L7x\J=2 
[29] PIMEO+1r21 
[30] CYCLE+CYCVECT[LT-1] 
(31] LAM+LAMSO2C 7-1) 
[32] CI+LAM NCTC 
[33] DFERVIAT+(CI-(LAM NCT(C-1E°8)))+#1F 8 
[34] MFWTOPRAPS2 
[35] +L7x\DIVERG=0 
[36] KILL«+1 
(37) C#CSTARCI ;J-1] CINCRNEWRAP2 ALP 
{38] KXILL<+<0 
[39] L7:CSTARCI3J1]+C 
C40] PTIMES<(( 121 )-TIMEO)+60 
met). ! 'sZOOP;! ';[ TIMES ;! 'sLITTRC IJ -11;'° ! 
C42] MULT+1 
C43] +Z14x1,JS#2 
C44] MULTH#(LITTROJI-LITTRO -1)) tC LI7TRC SJ -11-OLDCARRI133)) 
[45] +215 
C46] L14:4b15xr1is( (LIT TPL TS -11<50,4)+(( 7-1) =p L77TR)) 
(47) MULT+(LITTROSIJ-LITTRC I-11) e(LITPRO ST -1]-L7T7 TPCT - 21) 
C4e) 515 :C+C+eMULTx (C-CSTARCI 37-11) 
C49] £13 :+L11x%1(2+pLITTer )eJ+d+1 
[50] CSTAR(I311]+DFVECTCI]) 
(Si) '----- 22 22 - on nn nn oo on ee ee ee ' 
(52] 3+1%1(1+pDFPVEC™)#I+I+1 
(531 ™FMPLITT,R+0 ,LITTR 
[54] ARRAY+,CSTAR 
[55] APRAY*+TEMPLIT™TR ARRAY 
{56] ARRAY+(I, (7-1) )p ARRAY 
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MMIC ue one 


9 M+LAMV TRUNERR ERROP;7UNK3ELAM3N3K3E2R 
(1) AGLOBAL INPUT: LITTR,FRROR 
[2] TUNK+(K+1),(4+10) ,(LANMSQ2+0,5xLAMV x2) 
C3] Drs+b3xriLITTR( XJ >1 
cy *L11x .VLAMVEK1<6 
(5] oh4,(N+1),(CSTARCI sX+1]+ 1) 
C6] T11:+Z23*x 1.LAMV(X]20,1 


C7] Net 

C8] >L4 

[9] ba; te (CCP LAMVEX 1) #2) 42%. LAMVEX]) - -2 
(10) FLAM+LAMSO2(K)]x2,71828183 


(11) L2:EPR+(ELAM#N +1) #N44 
[12] M+N+1 

(13] +L2x\FERR>ERROR 

(14) LosMeM,N-1 

C15] +L1ixri(1+pLAVVy)#k+X+1 
. Vv 


VRP IB Q3sV3sM3A3sX323BP373003I3P3J3;PPECIS 
C1] aGLOBAL INPUT: XX(DF+ (DF+C#2)), PPFCISION( EQUALS 1+tHRROR ) 
[2] X*V/WXxXV*XXS0,5 
C3] J+(T+0),(R+0) 
C4] L1 :A+(Q=19)+(0-LO)xQ21Q 
C5] I+I+1 
C6] >Lux1(+/V)=0 
C7]. MefPf(@PRECISION:(P!A+P-1)x([/X)*P)4+0,6931471806 
C8] B+O 
C9] +L3x \A=1 
(10) BB+eliM)o, aM 
[11] BHM M)p((1M)-A) FM 
[12] B+«x/(BxBR)+1-BB 
[13] BHC H*«(X*P)xP!A+P-1)°o,*xRxPtPr yy 
C14] BeH++/BxXo,*xyM 
C15] L3:7+(0Q,Q0 )p(Q+1-100)P+0-190+L0 
C16] RB+(190)°,219Q 
(17] BR+«x/(HxRB)+1-BB 
[18] BR+((COYP+O-1)xX%P)°,xBR)x(1-X%)0, *Q-10Q 
(19] B+B++/BB 
[20] R+R+(V\BxT=1)+(V\1-B8)xI=2 
(21] L4:4+P 
[22] J+(V+~V) ,(O+H),(P+Q) 
(23) X+V/1-XxXx 
C24} +61%1(0<pX )ar=a1 
V 


UV CESTIMATE*+NEWR1 CESTIMATOR PRFEVIOUSAPR 3LAST2 37UNK 3I3DP°OLD sMUL™ 
(1) AGLOBA¥Y INPUTS+DFVECT ASSOCIATED WITH NEW LITTR | 
[2] JUNK+(Ret to PPEVIOUSARR ) ,(I+1),(CESTIMATE+ 10) 
C3] LAST2«(P, 2)+PREVIOUSARR 
ful DFOLD+, 1+PREVIOUSAPPL ; BU] 
[5] MULT#(NEWR1 -LAST2C1;3 21)8(LA AST20.1321-LAST20 1511) 
C6] 21:+b2%1124(+/(9PVECT=DFOLPLTI)) 
[7] CESTIMATE+CFSTIMATE, (MUL Tx (LAS T2C T+1321-LAST20 I[+1311]))+DAS720I+1327 
C8] £2:+L1%1(1+pDFOLD )zI+I+1 
Vv 
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SYNTAX: PTAIL*+DFVECT? NCT CVECT 


TRIS FUNCTION CALCULATES THE UPPER AFD LOWER TAIL PROBABILITY 
OF A NON-CENTRAL T DISTRIBUTED R,V, RELATIVE TO A DrSIRFD 
VALUE C(FOR FXAMPLE,PTAIL =PROBE|R,V,|>C]),9"°VECT IE A VECTOR 
(OR SCALAR) CONTAINING DFGREES OF FREEDOM OF INTEREST,CVECT IS A 
VECTOR(OR SCALAR) CONTAINING C VALUES OF INTEREST,FUNCTIONS ASSOC- 
TATED WITH(AND CALLFD AUTOMATICALLY) APF IB( INCOMPLETE RETA) AND 
TRUNERR( DETERMINES LEVEL TO WHICH THE INFIMITF SEPIFS ASSOCIATED 
WITH THE NCT CDF MUST BE SUMMED IN ORDER TO RFST?ICT THE TRUNCATION 
ERROR TO THE SPECIFIED VALUE).THF CPU TIME REQUIRED TO EXECUTE 
THIS FUNCTION FOR ONE VALUE IS GREATLY FFFFCTFD BY THE SPECIFIED 
VALUES OF DF,ERROF,C,AND AC(AS DEFINED BELOW),THF LARGER THE DF AND 
A,AND SMALLER THE ERROR THE GREATED THE CPU TIME,AS A POUGH IDEA- 
FOR ERROR+,0001,,01SAS,2,1<5DFS516; ,01<CPUS2 SFC, PER VALUE,FOR 
ERROR+,000001,,2<AS51,16SDFS50,2<CPUS120 SEC, PFR VALUE,TO COPY 
ONLY THF FUNCTIONS NECESSARY FOR NCT,COPY GROUP NONCENTRALT 


RESULTS WITH A TRUNCATION ERROR OF LESS THEN ,0001 FROM THF 

TRUE TAIL AREA) 

2,A+LOCATIOMN( DISTANCE) OF THE NCT FROM ZERO(NOTE: THIS DETERMINFS 
THE NONCENTRALITY PARAMFTER -LAM-), 

3,OPTION*+0 OR 1(0 WILL READ IN GLOBAL -A- AND CALCULATE -LAM- INTERNALLY 
»1 WILL READ IN -LAM”=- GLOBALLY AND CALCULATE -A- INTERNALLY FOR 
EACH DEGREF OF FRFEDOM 

4,DLAM*+NONCENTRALITY PARAMETER=Ax(DF+1)*,5(NEED ONLY BE SPECIFIED IF 
OPTION+1] 

5, HFADING+O OR 1(0 WILL SUSPRESS PRINTING OF A DESCRIPTIVE READING 
ABOVE THE OUTPUT PTAIL ARRAY),NOTE:SPACING FOR THE DFSCRIPTIVE 
HEADER IS DESIGNED FOR )DIGITS 4,5, OR 6, 


QUTPUT 
OUTPUT FROM THE FUNCTION WILL CONSIST OF AN APRAY( SCALAR) OF VALUFS 

WITH CVECT AS THE FIRST ROW,DFVECT AS THE FIRS” COLUMN, -A- AS THF 

SECOND COLUMN,-LAM- AS THE THIPD COLUMN, AND THF ASSOCIATED PTAIL IN 


THE REMAINDING COLUMNS AND ROWS 


V CYCLE+LAM TRUNEPP ERROR: TUNK;ELAM3N;ERR 
C1] JUNR+(CYCLE+((L LAM) *2)+3xLLAM) , (LAMSO2+0,5xLAM#2) 
[2] FLAM+LAMSO2%2, 71828183 
C3] Li:ERR+(ELAM!CYCLE+1 )*xCYCLE+1 
C4) CYCLE+CYCLE+1 
[5] +L1x1FRPR>ERROR 
(6)  CYCLE+CYCLE-1 

v 
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